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C D ' We developed a theoretical model able to give a common origin to the correlations 

, between the mass M, of supermassive black holes and the mass, velocity dispersion, 

kinetic energy and momentum parameter of the corresponding host galaxies. Our model 
ff^ ' is essentially based on the transformation of the angular momentum of the interstellar 

material, which falls into the black hole, into the angular momentum of the radiation 
emitted in this process. In this framework, we predict the existence of a relation of the 
form M, oc RcO^, which is confirmed by the experimental data and can be the starting 
point to understand the other popular scaling laws too. 
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H ' 1- Introduction 



At present, thanks to the improved angular resolution of modern telescopes, the ex- 
perimental evidence indicates that a variety of nearbygalaxies host a supermassive 
black hole (SMBH; M, > IO^Mq) at their centeJInSI. The subsequent discovery 
of a large number of scaling laws, in which the mass of SMBHs correlates with 
the properties of the host galaxies (bulges Jp^l^, demonstrates a link between the 
process of accretion of SMBHs and the formation and evolution of their galaxy. 

Even if many analytical and numerical models have been proposed at the same 
time in order to explain the observed scaling relationships (see for example Refs.fT2l- 
[T7| . the physical origin of these correlations, as well as the answer to the question 
"what is the most fundamental one?" , are still unclear and under debate^^'~'J^. 

Here we propose a simple theoretical model able to give a common origin for the 
correlations between the mass M, of SMBHs and the mass, velocity dispersion, ki- 
netic energy and momentum parameter of the corresponding bulges. Starting from 
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a principle of conservation of the angular momentum and usin g, a.s a first approx- 
imation, the Bondi-Hoyle-Lyttleton (BHL) theory of accretior|2QH22|^ found a 
fundamental equation of the form Af, cx RqCt^ , where Rc and a are the bulge 
effective radius and effective stellar velocity dispersion, respectively. From the pro- 
jections of this fundamental plane, using suitable correlations, we easily derive the 
other popular scaling laws. Despite the drastic hypotheses and hence the simplicity 
of the model, our results show an excellent agreement between its predictions and 
the experimental data, indicating that we have found a basic law for galaxies and 
their SMBHs. 



2. The model 

Marconi and Hunt in 2003 were the first to note that M, is significantly correlated 
both with cr, and with i?c- The conclusion of their study was that a combination 
of cr and Re is necessary to drive the correlations between Af, and other bulge 
propertieJ^. This topic was then deeply investigated through simulations of major 
galaxy mergers, which defined a fundamental plane, analogous to the fundamental 
plane of elliptical galaxies, of the form Af, c>c Rl^'^a^ or Af, c>c M^^'^a^, and Af, oc 
{Mi,a'^)°-'' , where Af* is the bulge stellar masJ^SHMI, These scaling laws are very 
similar to what was really found observationalljElHlHI, Following this path we build 
a theoretical model able to explain all the famous relations linking the mass of 
SMBHs with the properties of their bulges. 

Our model is mainly based on two hypotheses: the conservation of angular mo- 
mentum and a suitable velocity field of the gas. Then we need to make an approxi- 
mation to estimate the accretion radius of the black hole, which can be found either 
in a rough way or recurring to the BHL theory of accretion. 

Let us consider a black hole of mass Af, emitting radiation at rate L^, and 
accreting from a distribution p of gas, whose inward velocity is Vin- 

First, we assume that the angular momentum is conserved in such a way that 
a part (Afacc) of the total mass of the gas contained in the galactic bulge will 
be captured by the black hole, converting its angular momentum into the angular 
momentum of the perpendicularly emitted radiation (with an effective mass Afiad)- 
Of course, there are other mechanisms of conversion, transport or dissipation of 
angular momentum (viscosity, change in the spin of central black hole, etc.) that, in 
the past, have contributed with a different weight and importance to the accretion 
of black holes and probably are still acting now, but our drastic hypothesis focuses 
the attention only on the above described process. More complicated and detailed 
models should involve effects due to the other mechanisms that in our approximation 
are neglected. 

All orde r of magnitude estimate of the angular momentum for the galactic 
bulgJUHSO! leads to the conservation equation: 



Mace Rc Vrot = C Afjad Ra, 



(1) 
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where i?o is the effective radius of the bulge, Vrot is the mass-weighted mean rota- 
tional velocity of the gas, c is the speed of light and Ra the accretion radius of the 
Black hole. 

Then, we suppose that the velocity of the gas in the galactic bulge is related to 
the effective stellar velocity dispersion a in such a way that T^n = c (see Refs. [TT] 
and [31]) and Kot — Aa, where ^ is a constant!^. In the ideal case of the isothermal 
sphere, A is equal to ^/2. 

Now we must estimate the accretion radius and as a first approximation we can 
use the Bondi-Hoyle-Lyttleton theory " - : In their model the rate at which the 
gas is accreted onto the black hole is 

Mace = 7ri?| Finp, (2) 

where 

Rb = 2GM,/Vl (3) 

is the BHL radius and G is the gravitational constant. Of course the BHL theory is 
a realistic model in the case of radial accretion with low transverse velocity but it 
can be used as a good approximation even in presence of angular momentum (see 
the Appendix A) . Anyway, also without considering the BHL theory, the expression 
of i?B in Eq. (3), with Vln — c, roughly corresponds to the well known gravitational 
radius of influence of a black hole, which we can adopt as an estimate of Ra. 

Hence, we consider Ra — Rb (see the Appendix A), and substitute Eq. (3) 
in Eq. (1). Finally, recalling our hypothesis on the velocity field, we obtain the 
fundamental equation for supermassive black holes we were seeking: 

where e is the mass to energy conversion efficiency, which is set by the amount 
of rest mass energy of matter accreted onto the black hole that is extracted and 
radiated outward (e = Afrad/A/acc). 

The accreting black hole liberates energy at a rate 

eMaccC^ = ^^Edd. (5) 

It is commonly assumed that the accretion of matter onto a black hole releases 
energy at 10% efficiency so that a fixed value of e = 0.1 is usually adopted^Sl^ even 
if the range 0.001 < e < 0.1 has been recently investigated^^. 

This radiated luminosity is related to the Eddington limit, iEdd- In particular, 
ioT £ = 1 the central black hole is radiating at its Eddington limit: 

iEdd = ^ = . (6) 

CT KEdd 

Here mp is the proton mass, ctt the Thomson scattering cross-section of the electron, 
and KEdd = ct /"^p the opacity of the fully ionized hydrogen. Substituting Eq. ([2]) 
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and Eq. ([5]) in Eq. ([S]), and recalling our two hypotheses, we can determine also the 
gas density at the effective radius in a form: 



^Edd ^ Re ' 

similar to that found in Ref. [T31 While in other models the gas density is given 
among the hypotheses, in our approach it is a consequence of the theory. 

We have tested the effectiveness of our model on a sample of 58 nearby galaxies 
(z ~ 0). By using these experimental data and the corresponding errors, we report 
in Table 1 the best-fitting values for the slope m and the normalization h for the 
linear re lati ons used in this work. These values have been calculated by the routine 
FITEX^ for the relation y — b + mx, by minimizing the x^- The estimates of the 
Xied ~ X^/(58 — 2) and the Pearson linear correlation r for each relation are also 
shown. Errors have been calculated by using the formula reported in Appendix A 
of Ref. 22j The results of the fits for the M, — Mq and M, — Mqct^ relations are in 
agree ment with the ones obtained in a previous paper using three different samples 
of dataP^ and with the results obtained by other authors while the M, — a relation 
requires a more careful discussion (see Appendix B). The best-fitting line for the 



Table 1. Black hole-bulge correlations and fitting parameters for the 
considered galaxy sample. 



Relation 




m ± Am 




CO 


r 


M. - Roo-|oo 


8.11 ± 0.03 


1.00 ±0.02 


5.00 


0.39 


0.87 


M. - Mgo-2/c2 


4.56 ± 0.10 


0.87 ±0.02 


6.03 


0.40 


0.87 


M, — (7200 


8.21 ± 0.02 


5.83 ± 0.15 


6.22 


0.40 


0.87 


M, - Mgct/c 


0.73 ±0.19 


1.01 ± 0.03 


6.61 


0.42 


0.87 


M, - Mq 


8.74 ± 0.03 


1.21 ± 0.03 


7.71 


0.45 


0.85 


Rc — ly 200 


0.12 ± 0.01 
-1.77 ±0.06 


2.72 ± 0.09 
0.45 ± 0.01 


10.31 
3.26 


0.32 
0.16 


0.75 
0.92 



new relationship M, — i?oC2oo ^ log-log plane is: 

logio = (8.11 ± 0.03) + (1.00 ± 0.02) logio(i?caioo), (8) 

where (T200 is the velocity dispersion in 200 km sec~^ units, while is expressed 
in kpc. The linear relation ([8]) has a slope equal to the unity, which is exactly the 
value predicted by our model in the Eq. 

We remark also that this relation has the best and r among the relations 
in the upper part of Table 1 involving the black hole mass . We also report the 
intrinsic scatter eo, finding that the relation Af, — RcO^ has = 0.39, whereas the 
other ones have values of eg > 0.40. 

Furthermore, from the normalization we can calculate the effective efficiency 
coefficient, which turns out to be equal to i = 2£/A ^ 0.048 ± 0.003, a value close 
to 0.05 estimated in Ref. 112} and to 0.06 used in the case of Schwarzschild's metriJ^Sl 
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Of course if the coefficient e depends on Re or cr the relation (8) will still exist but 
its slope will be different from the unity. 

In Figure 1, we report the Af, — -Roflnn diagram in a log- log plot (we associated 

OA 

a particular symbol to each galaxy according to its morphology ). The best-fitting 
line is also shown. 

We also tested Eq. (|4]) by using an old sample of 37 objectJSl^ obtaining a slope 
equal to 0.90 ± 0.04, which is in good agreement with the slope of Eq. (8). Using 
a subset of 27 galaxies of the same old sample, Hopkins et al^^ obtained a fun- 
damental plane M. cx R"^ '^^cr^ °^^ and Graham^^ M, cx R'^ '^^a^-^^ . The difference 
with our result is due to a different fitting method and to the fact that these au- 
thors have considered three free parameters. With a different sample also AUer and 
RichstonJI^ studied the same three parameters fit obtaining M, cx R^ '^^a'^-^^ . 

The fact that we predict a relation of the kind M, cx Rcd^ and not the Hopkins 
result^^ Af, cx R^-^a'^ does not depend on an oversimplification of our model. 
To obtain the latter relation one must start from different hypotheses and only 
increasing the number of experimental data it will be possible to distinguish the 
better approach. 




Log,o[R<,o-2oo] 



Fig. 1. The Mm — Rofloo relation for the galaxies of the considered data set, where (T200 is the 
bulge velocity dispersion in units of 200 km s~^. The symbols represent elliptical galaxies (red 
ellipses), lenticular galaxies (green circles), barred lenticular galaxies (dark green circles), spiral 
galaxy (blue spirals), barred spiral galaxies (dark blue barred spirals), and dwarf elliptical galaxies 
(orange ellipses). The black line is the line of best fit for the sample of galaxies considered. 
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3. The origin of the other scaUng relations 

The dynamical masses of bulges, Afc, can be estimated by 

Mg = —Q—^ (9) 

where A: is a model dependent dimensionless constant. We stress that the isother- 
mal model (i.e. cr is a constant throughout the galaxy) is not a hypothesis of our 
framework. We make use of this approximation in order to test our model. As in 
Refs. |6] and |34j we use fc = 3 to compute the "isothermal masses" of the galaxy 
sample considerecP^ (in a more detailed model k is dependent on the Sersic index) . 
Thanks to the fundamental Eq. (|4]), coupled with Eq. in the following we de- 
rive all the most famous scaling relationships between the black hole mass and the 
parameters of the host galaxy. 

3.1. The M, — Mg (t/c relation 

Replacing Eq. (0) in Eq. (|H), we get: 

in optimum agreement with the corresponding relation in Table 1, where the value 
of the slope 1.01 is close to unity as we expected (the quantity Moa/c is also known 
as the momentum parameter-'^"'^). If we impose the exponent of the momentum to be 
equal to 1.00, then by refitting the data we obtain a normalization b — 0.82 ± 0.02, 
from which we derive e — 0.05. 

3.2. The M, — a- relation 

Substituting the experimental relation between Re and a (see Table 1) in Eq. (jl]), 
we obtain: 

M, = (a2oo)"' = 10'-'' {a,oo f'' , (H) 

ecG 

which is in agreement, inside the errors, with the corresponding law in Table 1: 

M, - 108-2i±o-02(^2oo)'-''^°-''- (12) 

3.3. The M, — Mg(t^/c^ relation 

Using Eq. Eq. Q can be written in terms of the kinetic energy of random 
motions: 

Replacing R^ in Eq. 1^ with the value R^ = 10"i-^'^(MGcrVc^)° '"'> taken from 
the experimental relation in Table 1, we obtain: 

M. = 104-^°(MGfTVc')°-'', (14) 
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where we have used the value of 0.048 for e. We point out the good match between 
the relation in Eq. (|14p and the corresponding one reported in Table 1: 

M. = 104-56±0.10(^^^^2/^2)0.87±0.02^ (^5) 

As shown in previous paperJ^^IMl^ this relation has the best compared to M, — a 
and M, — Mq laws. 

3.4. The M, — Mg relation 

Starting from the Eq. ([T0|) and expressing a in terms of M, (see Table 1), we obtain: 
in excellent agreement with the relation: 

M. = 10-4-57±0.34^^1.21±0.03^ (^7) 

4. Summary 

To sum up, we have shown that the model proposed in §2 works very well since, 
given a consistent set of data, it perfectly predicts the slope of a new relation (Eq. 
E]). Moreover, the most investigated scaling relationships can be easily obtained as 
projections of the plane identified by the fundamental Eq. Unfortunately, we 
cannot infer this correlation at redshift z >> 0, because we are limited by the small 
number of observable hosts. Other mechanisms of accretion, different from the one 
presented in this work, may have acted in the past, and hence other fundamental 
equations may have ruled the first Gigayears of the life of SMBHs. Future detection 
of new SMBHs, especially at higher redshift, and measurements of their masses will 
enable us to confirm the universality of our law or if it holds just for a restricted 
period of the cosmic time. At this stage it is early to predict the consequences of our 
equation in the context of the models about the co-evolution of galaxies and black 
holes. First of all we must check the validity of the approach using an enlarged 
sample of data. Recently Sani et al. collected a new interesting (but not larger) 
galaxy samplJ^Sl and we intend to present a complete analysis of their data in a 
forthcoming paper. A preliminary result is that a tight relation A/, — RqCt^ exists, 
but its slope oscillates from 0.78 to 1.2 depending on the number of pseudobulgecl 
considered in the fit. This aspect requires a deeper investigation before drawing any 
conclusions about the evolution of supermassive black holes and galaxies. 
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dynamics 
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Appendix A 

We underline that our model is not based on the BHL theory assumption. We can 
simply estimate the value of the accretion radius recurring to the concept of radius 
of influence of a black hole. Actually in our paper the BHL theory is equivalently 
used as another approach to have an estimate of the accretion radius even if it is 
well known that in presence of angular momentum the BHL approach is only a fair 
approximation. 

In their original paper, Hoyle and Lyttleton suppose that an element of volume 
of the gas cloud has an initial angular momentum, but it "Zoses this momentum 
through its constituent particles suffering collisions" at a radius Rb ■ That is, the 
transverse velocities of the particles, which reach the accretion line from opposite 
directions, annihilate reciprocally, whereas the radial component, if it is not bigger 
than the escape velocity (and this occurs at a distance d < Rb from the hole), 
makes possible that the particles were captured by the black hole. Then they show 
that the collisions occur with sufficient frequency to be effective in reducing the 
angular momentum. 

Anyway, it is possible to take into account the angular momentum, but the final 
result does not drastically change. In fact, in a generalization of BHL approach, 
the author in Ref. [22] considers the orbital velocity as the source of a lateral pres- 
sure on the accretion column. Having included this pressure term in the dynamical 
equations, he obtained an accretion radius Ra of the same order of magnitude of 
the BHL radius: 0.3 i?B < < 1.75 i?B. Both BHL and Horedt analyzed a case 
of pure accretion of the central object without considering the energy eventually 
radiated outward. 

Our framework is also different compared with the spherical Bondi accretiorP3, 
where a sound speed of the gas is considered, and a numerical correction factor ae ~ 
100 (see Ref. [T4| . which depends on the mass profile and gas equation of state, is 
inserted in Eq. ([2]). Some authors, by analyzing the role of angular momentum, refer 
to that special case of spherical Bondi accretion, finding that it fails systematically 
to reproduce their numerical simulations^. On the other hand, other authors, by 
using Chandra X-ray observations of several nearby elliptical galaxies, observed a 
tight correlation between the Bondi accretion rates (calculated from the observed 
gas temperature and density profiles as well as from the estimat ed b lack hole masses) 
and the power emerging from these systems in relativistic jetfP^J. They concluded 
that the Bondi formulae provide a reasonable description of the accretion process 
in these systems, despite the likely presence of angular momentum in the accreting 
gas. 
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We can study the effect of making the BHL theory a real assumption of our 
model. Let us analyze a different case (based on different hypotheses) from the one 
considered in section 2, assuming that the angular momentum is very small so the 
BHL theory is not only a fair approximation, but it works very well. Actually we can 
suppose that the photon is emitted by the black hole with a velocity c in the same 
direction of the arriving gas particle on the accretion line. It means that c forms an 
angle a with this line and with the radial component of gas velocity Vln and hence 
A = Vrot/Vin = tan(a). In this case the conservation of angular momentum can be 
written 

Afacc Rc Vy-ot = Mace V^n tan(a) = c Mrad i?A sin(Q;). (18) 

If we suppose that the accretion occurs mainly for particles with low angular mo- 
mentum (A << 1), we can simplify the above equation because for a ~ we have 
tan(Q!) ~ sin(Q;). In this particular case, our fundamental Eq. (4) becomes 

and the effect to consider the BHL theory as a real hypothesis of the model reduces 
to the disappearing of the factor A from the fundamental Eq. (4), and all the 
calculations of the first part of the paper are still valid considering e = 2e. We 
reported again the data of the Hu samplJ^H in Figure 2 together with the lines 
representing the values of the efficiency coefficient e predicted by Eq. (19). So, 
given the SMBH mass and the effective radius and dispersion velocity of the host 
galaxy, the equation (19) allows a quick estimate of the efficiency of a black hole. 
In this way the resulting diagram (Fig. 2) can be used to classify the black holes 
in terms of their efficiency provided that the host galaxies satisfy the hypotheses of 
the model. For example, in our sample all the galaxies have a predicted e < 0.25, 
but the prediction could be not valid, as expected, for some pseudobulges for which 
the approximation of the model ^ << 1 does not work. 



Appendix B 

The slopes of the M, — Mq and M, — Mqo^ relations remain very stable changing 
the sample of data or the fitting method. In Ref 28 we have used one fitting method 
and three different samples and the resulting slopes were (1.18, 1.22, 1.27) for the 
M, - Mq relation and (0.83, 0.86, 0.91) for the M, - Mqu"^ relation. Viceversa, 
in Ref. 27 we have used only one sample and three different fitting procedures 
and the slopes are still stable with values (1.15, 0.98, 1.07) and (0.80, 0.74, 0.78) 
respectively. These values are in agreement with the results of other authors. For 
instance, Hopkins et al.^2Slfin(j i_05 for M, — Mq and 0.71 for M, — MQa"^, whereas 
Soker et al.lJJJfind 1.07 and 0.74 respectively. 

The M, — a relation has a different behavior since its slope strongly depends on 
the fitting method (5.06, 4.46, 4.25) in Ref [271 and on different samples (5.26, 4.99, 
5.83) in Ref. 1281 In this case Hopkins et al.'^nd 3.96, SokeJIIl4 .18, Ferrarese and 
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ForcPni4.86 and 5.1 and HiJSl a slope that changes from 4.01 to 5.62 depending on 
the considered subsample of data. 

The oscillation of the slope between 3.9 and 5.9 could be due in part to the 
fitting method (FITEXY finds the first minimum of and gives a slope; if one 
introduces the intri nsic scatter and refits until the reduced = 1, a smaller slope 
is generally obtained^^l) and in part to the presence of pseudobulges in the sample. 
In their last paper Sani et. al.l^find the same value (4.00) for the slope using three 
different methods, but selecting only the classical bulges. By introducing also the 
pseudobulges in their sample, the slope increases and the result becomes no longer 
stable. 
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